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The prospect of coherent dynamics and excitonic delocalization across several light-harvesting structures in 
photosynthetic membranes is of considerable interest, but challenging to explore experimentally. Here we 
demonstrate theoretically that the excitonic delocalization across extended domains involving several light¬ 
harvesting complexes can lead to unambiguous signatures in the optical response, specifically, linear absorption 
spectra. We characterize, under experimentally established conditions of molecular assembly and protein- 
induced inhomogeneities, the optical absorption in these arrays from polarized and unpolarized excitation, 
and demonstrate that it can be used as a diagnostic tool to determine the resonance coupling between iso- 
energetic light-harvesting structures. The knowledge of these couplings would then provide further insight 
into the dynamical properties of transfer, such as facilitating the accurate determination of Forster rates. 

Keywords: Photosynthesis, Excitons, Quantum Delocalization, Purple Bacteria, Spectroscopy 


I. INTRODUCTION 

Nature has evolved a variety of photosynthetic archi¬ 
tectures. A detailed, quantitative understanding of the 
principles that underly their function could assist the de¬ 
sign of future energy conversion devices. A wealth of 
careful structural and spectral studies 1 ^ complemented 
by first principle calculations^" 13 - have contributed to our 
current understanding of exciton transfer dynamics in 
light-harvesting (LH) antenna and reaction center (RC) 
pigment-protein complexes. Recently, interest in this 
topic has intensified due to observations of persistent os¬ 
cillatory features in non-linear optical experiment^®®} 
which have been reproduced in various LH structures. 
This has motivated work which reevaluates the nature 
of the interaction between excitonic dynamics and vibra¬ 
tional motiorP®®l. 

Recent work has examined excitonic delocalization 
restr icted to closely-packed pigments in single LH 
unit^ 22 23 . This approximation is valid for studies of 
inter-complex energy transfer and fluorescence, which 
take place on time scales (~ ps-ns) much longer than the 
dephasing (~ 100 fs), which re duce s the excitons to lo¬ 
calized (single-unit) domain^- 10 * 24 * 25 * However, as we will 
show, for processes that are characterized through ob¬ 
servables with a faster built-in time-scale, such as absorp¬ 
tion, which depends on the dipole-dipole correlation func¬ 
tion (DDCF) oscillating at optical frequencies, extended 
delocalization across LH complexes must be taken into 
account in an accurate description^. Here, we propose 
and characterize theoretically an experimental scheme, 
based on simple linear absorption measurements, which 
quantifies the resonance coupling among LH complexes 
in purple bacteria. Based on firm theoretical analysis, we 
show how this information facilitates the determination 
of incoherent excitonic transfer rates in systems where 
the donor’s fluorescence and acceptor’s absorption over¬ 


lap. 

II. MATERIALS AND METHODS 

Purple bacteria. Photosynthetic membranes of pur¬ 
ple bacteria are composed, in general, of two-dimensional 
arrays of LH1 and LH2 complexes, which are responsible 
for the absorption of light and the subsequent transport 
of excitation energy to RC pigments, where charge sep¬ 
aration dri ves met abolism under photosynthetic growth 
conditionIn such membranes, excitonic en¬ 
ergy transfer among LH1 and LH2 units has been char¬ 
acterised in terms of incoherent (thermal) hopping 31 . 
However, for a description of absorption, which has a 
much faster built-in time-scale than inter-complex energy 
transfer and dephasing, excitonic delocalization across 
several harvesting complexes may become relevant. In 
this work, we analyse excitonic delocalisation across 
many LH complexes. We will show that this delocal¬ 
ization leads to a redistribution of absorption intensity 
which is experimentally accessible by exploiting the sym¬ 
metry of LH1 and LH2 complexes. 

For concreteness, we develop our analysis using the 
LH1 of R. rubrum, whose assumed homology of LH1 to 
LH2 harvesting units, permits a straightforward general¬ 
ization to membranes composed of either complex. The 
available X-ray structure of the LH1 complex from Tch. 
tepidurM 2 and vanishing fluorescence anisotropy of LHls 
in R. rubruiA 33 provide support for a closed ring struc¬ 
ture with dimeric repeating units of 2N = 32 bacteri- 
ochlorophyll (BChl) pigments in a Ci6-fold symmetry, as 
depicted schematically in Figure [TJ a), although not all 
LH1 complexes are closed or circular®. LH1 complexes 
from R. rubrum naturally aggregate via protein domain- 
mediated interactions into arrays with a tetragonal pack¬ 
ing and can be grown without the constituent RCip] Due 








2 




FIG. 1. (a) The Q y transition dipole moments (red arrows) of the BChl pigments and the protein scaffold (light red) of the 
LH1 complex in R. rubrum, are pictured schematically® 1 . The rate of incoherent energy transfer between excitons k and k! on 
different LH complexes, '•/k.k ,', is determined by the resonance coupling 14, as explained in the main text, (b) The oscillating 
dipole-dipole correlation function (DDCF) contains all signatures of the resonance coupling and decays due to decoherence, 
which leads to homogeneous spectral broadening. The difference between the x- and y-polarized correlation functions (c) 
contains the signatures of coupling between rings, 14 ,fc', and can be obtained experimentally by measuring (d) the linear 
dichroism ( LD ), which is the difference between polarized absorption along these orthogonal axes. Importantly, the effects of 
local dephasing on the LD can be reproduced by both the DDCF, Eq. 2, or a simple dressing of excitonic states according to 
Eq. |3] 


to symmetry, the complex develops a single bright band 
at 880 nm, which, in the case of R. rubrum , has 2.4 times 
the dipole strength from what is expected from the ad¬ 
dition of 32 dipoles Such hyperchromism is ac¬ 

counted for by increasing the magnitude of the induced 
dipoles Id! | —1 |\/2.4 x d\ = \/2Ax6 .3 Debye = 9.8 Debye. 
In R. sphaeroides , LH1 complexes tend to dimerize 2 " 4 . 

Exciton formalism. Excitonic properties of mem¬ 
branes subject to low-intensity illumination can be ob¬ 
tained from the electronic Hamiltonian in the single ex¬ 
citation subspace: 

Jf = E e m | m) (m\ + E J mn \n) (m\ (1) 

m m,n 

where e m is the excitation energy of pigment m and 
Jnm is the coupling, through the Coulomb exchange 
mechanism, between the Q y induced transition dipoles, 
di, of the electronic excited states | n) and | m) on pig¬ 
ments n and rrP. Details of the coherent interaction 
lead to specific eigen-frequencies uj a from eigenstates 
| a) = Yln c n\ n )i ln general, delocalized over the com¬ 
plete set of pigments accounted for (see Supplementary 
Information for Hamiltonian parameters). In general, the 
resonance coupling between different rings, labeled Vk,k' 
in Fig^a), is the origin of EET, which is well-described 
as an incoherent process between different rings (occur¬ 
ring at a rate 7 ^*/ oc V),././ | 2 )P^. The key prediction of 
this article, illustrated in Fig0 is that this resonance 


coupling introduces experimentally-accessible signatures 
in the linear absorption which evidence excitonic delocal¬ 
ization across different rings and can be used to obtain 
estimates of the EET rates r )k,k>- 

The absorption spectrum Ag(ui) along an excitation 
polarization direction E is obtained from the Laplace- 
Fourier transform of the DDCF, (( D(t ) • E)(D( 0) • E)f^, 


The propagator contains all excitonic couplings and 
dispersive processes of relevance for excitonic dynamics. 
Within the DDCF description, the evolution of the dipole 
moment operator D(t) = e Ct D{ 0) = e Ct JT di \m) (0| 
leads to oscillations at optical frequencies oc e lulat with a 
period of about ss 2-3 femtoseconds, due to the Hamil¬ 
tonian part of the propagator. These coherent signa¬ 
tures are degraded by environmental decoherence, which 
in the worst-case scenario (in terms of processes that de¬ 
grade the excitonic delocalisation) may lead to indepen¬ 
dent fluctuations of pigment energies uj t within the ex¬ 
citation’s lifetime, local pure depha sing, occurring on a 
time-scale of 74 1 ~ 50 — 100 f d 4 ° l 41 ! 

The DDCF evolves for many tens of optical cycles, as 
seen in Fig. 1(b), before dephasing sets in, which can 
be sufficient to imprint in the absorption spectra fea- 
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tures of excitonic coupling, i.e., excitonic delocalization 
among several LH units. Such features are difficult or 
impossible to observe in processes with a longer built-in 
timescale, such as EET. Our interest is to connect the op¬ 
tical signatures from this DDCF function with excitonic 
delocalization across multiple rings. In order to be con¬ 
servative, we initially consider the worst-case scenario for 
degradation of such extended excitonic delocalization, in¬ 
troduced by a local dephasing superoperator of the form 

£dP = E m - p). with = \ m > m - 1 °) < 0 | 

and the density operator p, consistent with random tem¬ 
poral fluctuations of the energy gap between levels | m) 
and |0). When derived from a microscopic approach, 
the considered dephasing model describes the action of 
a fully Markovian environment in the infinite tempera¬ 
ture limit. The vectorial nature of the dipole moment D 
permits discrimination of absorption along the axis that 
joins two coherently-coupled rings and perpendicular to 
it, based upon the calculation of the respective DDCF 
of Fig. 1(b). The subtraction of these two signals, pre¬ 
sented in Fig.1(c), shows that their DDCF differ, a fact 
made even more conspicuous by the subtraction of ab¬ 
sorption spectra (through the Laplace-Fourier transform 
from Eq. [2| corresponding to these orthogonal directions, 
namely the linear dichroism ( LD ). This result should be 
contrasted with the vanishing LD of single circular rings, 
or equivalently, assemblies of rings where the resonance 
coupling between rings plays no role, as discussed below. 
As these calcuations show, strong local dephasing with a 
rate of 7d=l/50 fs is not enough to smear out the finite 
LD for the coupled rings, and displays the robustness 
of optical measures for understanding the extent of ex¬ 
citonic delocalization in photosynthetic membranes. Re¬ 
assuringly, the result obtained from the local dephasing 
model is virtually indistinguishable from the usual proce¬ 
dure illustrated by Eq. [3j in which excitonic energies uj a 
are used as the center of a dressing function /(w — u a ), 
whose weight in the complete absorption spectrum is de¬ 
termined by the exciton dipole strength \D a ■ E | 2 . This 
agreement justifies the reduction of the complexity of the 
DDCF calculation to a simple diagonalization of the full 
Hamiltonian and the use of Lorentzian line-shape func¬ 
tions f(u>) = 7 2 /(w 2 + 7 ^), in order to establish the in¬ 
fluence of the electronic coherent coupling in polarized 
absorption spectra. Moreover, it emphasizes the neces¬ 
sity of the excitonic dipole strength |D a | 2 , and therefore 
the relevance/necessity of excitonic states |a), for the de¬ 
scription/understanding of the the absorption process in 
ring assemblies. 

III. RESULTS AND DISCUSSION 

Excitons | a) are affected by their interactions with 
the protein environment, which dynamically degrade the 
electronic coherence within the excitation lifetime (ho¬ 
mogeneous broadening) captured by either Cd or the 
dressing procedure, and quasi-static fluctuations (inho¬ 


mogeneous broadening) of the excitation energy and the 
Coulomb coupling. An interplay between both types of 
broadening leads to the observed width of the absorp¬ 
tion spectrum, which in the case of LH1 is approximately 
r « 465 cm -1 . If the homogeneous line shape function 
is consistent with a dephasing lifetime of 7 ^ = 1/50 fs, 
a standard deviation of a Gaussian pigment site-energy 
variation of 325 cm ” 1 is required to recover the full width 
of the LH1 absorption (further details in the Supplemen¬ 
tary Information). This la tter value is in agreement with 
previous estimation d 24 ! 42 * ^. 

Extended delocalization: Absorption and LD 
as a witness of delocalization in ring assemblies. 
The relevance of the excitonic dipole strength \D a \ 2 for 
the description of the optical response of ring assemblies 
motivates our study of the excitons |a) in multiple ring 
arrays. 

In order to estimate the influence of inhomogeneities 
on the e xtent o f excitons, the usual inverse participa¬ 
tion raticP -HUD is generalised to quantify delocalisation 
across rings R of a given exciton a, namely the ring par¬ 
ticipation ratio IPR% = 1/Em=i(E^ |c“| 2 ) 2 - The 
IPRr ranges from 1 to R, and IPRr > 1 unequivocally 
represents excitonic delocalisation over more than a sin¬ 
gle ring. An estimation of the size of excitons in each 
of the three configurations presented in Fig j2ja) shows 
that absoprtion in these arrays generates excitons that 
in general extend over domains larger than a single ring. 
This figure also shows that excitons have similar sizes for 
a fixed set of coherently coupled rings. In more detail, 
small differences emerge from the specific connectivity 
of each configuration, setting the square, T-like and lin¬ 
ear configurations with exciton sizes in descending order. 
Importantly, the histogram of IPRr of Figj2] is made 
upon the two most optically active states from every noise 
realisation, and therefore represents the accesible states 
through optical excitation. 

For uncoupled LH1 complexes, the circular symmetry 
leads to an optical response predominantly determined 
by two orthogon ally-polariz ed degenerate exciton states 
and vanishing GUSD However, the coherent exci¬ 

tonic interaction bettween rings may lift such a symmetry 
and result in shifts in the unpolarised absoprtion spec¬ 
trum and/or to a finite LD. The key theoretical predic¬ 
tions of this article are summarized in Figures |2^b)-(e). 
Namely, a shift in the unpolarised absorption spectra will 
generally be observed in assemblies of many rings, while 
a non-zero LD will be measured in arrays with a large 
aspect ratio and small width, i.e. linear-like arrays. 

A close inspection of the absorption spectrum 
(Fig j2jb) ) of the configurations from the inset of Fig|2ja), 
shows very similar spectra for the T-like and linear config¬ 
urations, while the square assembly does present a larger 
shift to the red than the former geometries, in compar¬ 
ison to the spectrum from a single ring. In practice, 
shifts that account for fractions of nanometers can be 
difficult to discriminate and a complementary measure, 
namely the LD, can be of help to understand whether 
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FIG. 2. Extended delocalization across multi-ring arrays in different multiple-ring configurations. Each four-ring configuration 
considered in (a) exhibits a similar, although not identical, extent of delocalization, measured by the ring inverse participation 
ratio IPRr. Plotted is the probability distribution, dp(IPRn), of the IPRr of the exciton state with greatest dipole strength, 
averaged over realisations of static disorder. The average delocalization length, (IPRr), increases slightly with the number of 
neighbors, i.e. (IPRr) = {1.24,1.27,1.33} for the line, T-, and square configurations, respectively. However, the delocalization 
is witnessed by linear dichroism (LD) only when the assembly exhibits an asymmetry, like in the T- and line assemblies, as 
shown in (c). (d) Red-shifted absorption and (e) non-zero linear-dichroism (LD) for linear arrays (Q= 1, calculated for different 
values of R = 1,..., 5) witness extended delocalization over multiple LH1 complexes. Results from 5 x 10 4 realisations of 
site-energy variations with standard deviation 325 cm -1 and jd = 1/50 fs _1 . 


the shifts observed are compatible with excitons formed 
across several rings. Even though all these structures 
present a finite shift with respect to the single ring spec¬ 
trum, notice however, that not all multiple-ring arrange¬ 
ments present a finite LD. For instance the symmet¬ 
ric square assembly- which presents the greatest IPRr - 
has vanishing LD, while the linear array shows a LD 
with a contrast of ~ 4 — 5% of the absorption maximum. 
From these results it is apparent that, although each as¬ 
sembly presents similar delocalization and shifts in the 
unpolarised absorption, asymmetric assemblies provide a 
suitable scenario for the experimental observation of exci- 
tonic delocalization using polarized spectroscopy. When 
the focus is placed in the linear-like arrays, it is observed 
that unpolarised absorption spectra shows greater shifts 
in Figj2^d) with increase of the number of rings in the 
linear array. Interestingly, the amplitude of the LD also 


increases with the linear array size, from ss 2.5% for 2 
rings, rising up until it seems to saturate to « 5% for 5 
rings. 

It must be highlighted that just due to homogeneous 
broadening, the LD amplitude for 2 rings is about 5% 
while if the inhomogeneities are included, it reduces to 
about 2.5% of the total spectra. We have tested that 
this latter value is robust to different compositions of ho¬ 
mogeneous and inhomogeneous broadening restricted to 
the same full broadening of the optical transition. The 
inhomogeneities in a sample may not be restricted to 
environmentally-induced inhomogeneities in the pigment 
energies and can indeed be caused by structural pertur¬ 
bations, such as variations in the inter-complex distance 
or geometrical perturbations that produce ellipticity. A 
thorough study of both effects is presented in the Sup¬ 
plementary Information, from where it can be concluded 
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that the specific details of the broadening mechanisms do 
not play an important role for the determination of the 
LD amplitude. 


These full numerical simulations of the polarised and 
unpolarised absorption spectrum highlight 1 ) in general, 
(small) shifts in the absorption spectra are expected in 
arrays of rings; 2) a finite LD , in particular for linear-like 
assemblies, will be observed due to the coherent coupling 
of the rings, and 3) the contrast of the LD depends in¬ 
versely on the total broadening of the optical transition 
and not on the specific details of the broadening mecha¬ 
nism. 


Analytical expressions. In order to under¬ 
stand which quantities determine both the shifts in 
the unpolarized spectrum and the LD , an analyti¬ 
cal model is desirable. In the absence of inhomo¬ 
geneities, the optical response of circular aggregates, 
resembling LH1 or LH2 rings, is determined by the 
\k = ±1) = EJ- 1 )” ex P (2nik[n/2\/N) | n) states, 

where |_ri/ 2 J is the largest integer smaller than n/ 2 . 
These states carry all the transition dipole moment, 

namely D k -±i ss \d\^J~^{x ± y), where x and y denote 
orthonormal axes to be identified with those shown in 
Figure pi a). The inter-ring resonance coupling Vk,k> = 
(k,e | Jr\k',e') between excitons k and k' on rings cen¬ 
tered at lattice sites e and e', calculated by coupling 
all pigments on both rings according to Eq. |T| is 
much smaller than the energy difference between single¬ 
ring excitons, namely | (k, e \ Jf’\k, e) — (k', e'\ \k\ e') \. 

Such condition results in a minor mixing between states 
\k\ ^ \k'\, results in an optical response of assemblies 
mainly determined by the coupling between the bright 


| k = ± 1 , e) and 


k' = ±l,e' 


states from neighbouring 


rings. In particular for this manifold, it can be shown 
that Vqi = V-i-i < 0 , Vi-i = V'* lil and the ar¬ 
gument arg(Vr,_i) = £ for e — e' = Lax becomes 
arg(y 1 _ 1 ) = £ + 7 r for e — e‘ ' = Lay , where a is the 
lattice spacing. 


In what follows, we consider nearest-neighbor cou¬ 
plings on the rectangular lattice of Figure [3ja), as in 
Figure [ 2 ] regarding square and linear configurations, as it 
exhibits the optical signatures of extended delocalization 
while allowing simple analytical expressions. In partic¬ 
ular, Vqi = -1.8 cm -1 and |Fl,_i| = 4.4 cm -1 . Qual¬ 
itatively similar results are obtained for couplings be¬ 
yond nearest-neighbours and triangular lattices, which 
we consider explicitly in the Supporting Information. 
Rings are centered at lattice sites e = a(qx L ry), where 
q = £ Z, r = 1 £ Z and a = 120 A® 

single-ring excitons are labeled |±l,y,r) = |±l,e). Ac¬ 
cordingly, the extended excitonic wavefunctions over the 


rectangular array read as 



(4) 

where |±) = ^ (|1) ± |—1)) and k q = 1, ...,Q and k r = 
are Fourier quantum numbers corresponding to 
the x and y directions, respectively. The optical response 
can be approximated according to equation [3] 

f\D + , ka , kr - E\*\ S\k qi Q)S 2 {k r ,R) fcos 2 <A ^ 
\\D_, kqikr ■ E\ 2 ) Q + l RL 1 Vsin 2 W l j 

where <f> is the angle between the polarization of the field 
E and the axis x. The absorption strength is determined 

by S(k,W) = sin (^) sin ^ 2 (WTI)) cs c ^ 2 (WTT)) > dis¬ 
tributed over states j±) along orthogonal polarizations. 
Given that S(k, W) ~ ^ sin 2 (^) already for W > 3 — 4 
rings, the dipole strength becomes small for higher k > 1 
states, concentrating dipole moment in the k q , k r = 1 
states in the rectangular configuration. The energies of 
these states are shifted by 

£ ± = 2y,, 1 (cos(^ T ) + cos( i ^)) 

±2|Vi,_,| (cos (gE) - ™ ( jAi)) 

= SE± A E. (6) 

Given that the absorption is composed of these two 
shifted optical transitions, a general shift SE = ( E + L 
E -)/2 oc V -[,1 of the (unpolarised) absorption maximum, 
with respect to the single-ring case, is expected for the co¬ 
herently coupled array of rings. Such a shift from the as¬ 
semblies studied in Figpla)-(b), accounts for SE = 2V 1:1 
in the square (Q = 2, R = 2) and SE = 2V\ t \ cos( 7 r/ 5 ) in 
the linear (Q = 4, R = 1) assemblies. Hence, the above 
expression explains why the square geometry shows a 
greater shift than the linear assembly, with respect to the 
single-ring case, in the unpolarized absorption spectrum. 
Last but not least, note that the magnitude of V 1.1 de¬ 
pends also on the specific details of the structure (further 
details in the Supplementary Information). Hence, spec¬ 
tral shifts can be used as a tool to understand specific 
details that have been not possible to resolve by other 
means. 

The states |+) and |—) are orthogonally polarized; po¬ 
larisation along the x-axis results in absorption spectra 
peaking at energy E + corresponding to D + oc x, while 
polarisation along the y-axis results in absorption spec¬ 
tra peaking at energy £L corresponding to the D_ oc y 
state. The subtraction of these two spectra recorded from 
perpendicular polarisations, the LD , is finite whenever 
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(a) 




FIG. 3. (a) LH complexes in purple bacteria naturally arrange themselves in well-ordered 2D arrays. The spectral consequences 
of extended delocalization in such an arrangement can be understood by examining a rectangular Q x R array. The absorption 
anisotropy is strongest for linear-like arrays ( Q/R 1), represented schematically by the color gradient, (b) The circular 
symmetry of single LH complexes concentrates the optical absorption strength in degenerate states with orthogonal polarization 
(blue bounding boxes). The circular symmetry of the combined system is broken upon interaction (orange bounding boxes), 
leading to a measurable splitting of polarized states. 


A E = (E + — EL)/2 does not vanish, which occurs, ac¬ 
cording to Eq. [Gj strictly due to the resonance coupling 
V k)k > which breaks the circular symmetry of the assem¬ 
bly. Explicitly, based on Eq. [6] the effect is greatest for 
a linear chain (Q=lorf? = l) and a saturation for 
both the red shift, SE , and splitting, A E, is expected, as 
cos gq^j- « 1 already for Q « 5 — 6 rings. 

Experimental determination of Forster rates. 
Besides controlling the linear optical response, the 
resonance coupling dictates the rate at which, on a 
longer time scal e, exc itations migrate incoherently be¬ 
tween complexe j — I— ( The transfer among different 
bands which can be resolved spectrally permits ex¬ 
perimental schemes that make use of the different 
spectral components in order to separate dynamical 
contributions^iEMsni j n the case of isoenergetic transfer 
steps, like LH1—>LH1, this is not possible and it becomes 


hard to identify the nature of individual dynamical com¬ 
ponents. Excitation annihilation was used to determine 
the transfer rate between LH2 rings in R. sphaeroide 
and techniques based on transient depolarization, such 
as anisotropy absorption recovery, have been valuable to 
recognize multiple dynamical contributions at low tem¬ 
peratures in the core complex of R. rubrurr |5®Hs31, but have 
failed to unequivocally identify components arising from 
intra-ring relaxation or from inter-ring energy transfer. 
The accurate determination of these rates is crucial for 
understanding the efficiency of photosynthesis. 

We show here that the possibility to obtain the res¬ 
onance coupling between isoenergetic species of LH1 or 
LH2 complexes- through the absorption measurements 
illustrated in Figure [3] circumvents the ambiguity in 
these isoenergetic landscapes and opens up a promising 
experimental scheme to quantify their mutual transfer 
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FIG. 4. Absorption measurements can provide direct information on the magnitude of transfer rates. In (a) is presented the 
transfer time between LH1 complexes as a function of the LD contrast and/or shift in the unpolarized absorption (which depend 
on the resonance coupling 14,*/, itself a function of the dielectric constant). An experimental ensemble measurement of the 
LD signal mitigates the effects of structural disorder and allows a direct absorption readout. Such an experiment is shown in 
(b), where “blue” rings, mutated to maximally absorb at less than 880 nrd 36 * 51 ( are grown in a membrane with “red” rings, 
absorbing at 880 nm and linked in pairs. The aggregation of blue-shifted LHls leads to their vanishing LD, while, if the axes 
of the red pairs align, their non-zero LD signal can be measured. 


rates. 

Two facts of major importance are that the shift of the 
unpolarised absorption is SE = Vi,i cos(7t/(Q + 1)) and 
that the splitting between bright tranitions is 2AE = 
4| Vi,_i| cos( 7 t/(Q + 1)) for linear geometries. Accord¬ 
ingly, measurement of the shift in the unpolarized ab¬ 
sorption spectrum determines Vj,i while, as will be shown 
shortly, the determination of V\ _i through an LD mea¬ 
surement is possible. A good quantitative agreement be¬ 
tween the full calculated spectrum and the subtraction of 
two Lorentzian lineshapes with full width T ss 465 cm ' 1 
(in agreement with experiment) separated by 14 ,_i per¬ 
mits to extract the relation 


LDmax — C\ 



(7) 


for the LD amplitude, LD max , with a constant C\ = Op 
in the case of two rings, which accounts for the geomet¬ 
rical details of the array. Accordingly, it is possible to 
determine through the unpolarized spectrum and the LD 
amplitude, the values of V \,-| and V/,-i independently. 

Generalized Forster theory is often used to calculate 
the rate of incoherent transfer "fk^k' from a donor exciton 


k to acceptor exciton k' 


Ik^fk’ = -r~|14,fe '| 2 Zk Ik,k', (8) 

where Vk,k' is, as before, the coherent excitonic coupling 
among single-ring states, Zk is the thermal population 
(Boltzmann factor) and Ik,k ' = f 0 Ek(e)Ak'(e)de de¬ 
notes the spectral overlap integral of the donor fluores¬ 
cence from exciton k and ac cepto r absorption of exci¬ 
ton k ', F and A, respectively^®!. The total transfer 
rate from donor to acceptor is 7 = Efc k’ 7 fc->fc' • Al¬ 
though we have taken a large inhomogeneous broaden¬ 
ing in our calculations, there is experimental evidence of 
a predominance of homogeneous broadening of LH1 at 
room temperatur^^®!. Under the assumption that the 
emission is also spectrally homogeneous, the total rate is 
7 ~ Efc.fc'=±i 7 k^k 1 , and the Forster rate can be calcu¬ 
lated as 

O 7 r 

7 = — (2|U 1 , 1 | 2 + 2|V 1 ,_ 1 | 2 )ZD (9) 

where / = J 0 °° F(e)A(e)de is the full spectral overlap 
and Zi is the thermal population of the \k = ± 1 ) states, 
which depends on the single-ring model. 











As shown above, Vi,i = -1.8 cm -1 and |Vi,_i| = 4.4 
cm -1 are quantified by shifts in absorption and non-zero 
LD contrast, respectively, which allow a calculation of the 
Forster rate with a minimum of model input. The dipole 
moments \d\ and the dielectric contant k = 2 predict an 
inter-ring transfer time of 8 ps, which is well within the 
bounds (~lps - 20ps) of estimat es from current experi¬ 
ments and theoretical calculations- flSSliHI Determination 
of the inter-ring coupling strength by measurement of the 
LD thus provides a new method to determine Forster 
transfer rates. This orocedure would then provide fur¬ 
ther insight into the values of the dielectric constant and 
the BChl dipole moment in vivo , which are parameters 
usually estimated indirectly. Figure [4] (a) plots the trans¬ 
fer time, LD contrast and absorption shift as a function 
of these quantities according to Eq. [9] 

Several elements have to be brought together in an ex¬ 
periment which confirms excitonic delocalization across 
extended domains. A speculative setup is illustrated 
in F igure [f] (b) based on mutagenesis of blue-shifted 
LH1^26|5U if a link between LHls by means of me¬ 
chanical bridges or affinity domains is accomplished, 
the formation of linear-like assemblies will be possible. 
Additionally, these assemblies must be macroscopically 
aligned in order to set the polarisation excitation direc¬ 
tions. The blue-shifted LHls will serve then to build 
the membrane scaffold to accomodate the linear arrays 
of non-mutated LHls and fix their orientations. Fluo¬ 
rescence yield detection upon polarised excitation per¬ 
formed on single assemblies® is also a promising tech¬ 
nique to confirm/disregard the hypothesis of long-range 
excitonic delocalisation in natural LH aggregates. 


IV. CONCLUSIONS 

In summary, we have modeled and characterized, 
through analytical expressions and numerical calcula¬ 
tions, the optical response due to resonance coupling 
between multiple rings in photosynthetic membranes of 
purple bacteria. This investigation highlights the impor¬ 
tance of a description, at room temperature, of light ab¬ 
sorption beyond the standard assumption of excitons re¬ 
stricted to individual rings. The dipole moment redistri¬ 
bution that emerges in 2D arrays of LH1 rings, due to ex¬ 
tended excitons, leads us to propose an experimental pro¬ 
cedure, based on polarized and unpolarized absorption 
spectra of small linear-like arrays of rings, which quan¬ 
tifies inter-ring resonance coupling. We show that this 
experimental procedure allows an alternative, indirect 
measurement of their incoherent Forster transfer rate, 
which carries additional information on currently poorly- 
characterized parameters, like the dielectric constant or 
the in vivo BChls dipole moment. Such a measurement 
could be accomplished through absorption of an ensemble 
of oriented linear-like assemblies or through fluorescence 
yield detection of single assemblies which have been ex¬ 
cited by polarized fields. 
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Appendix A: Hamiltonian of LH1 complex 



FIG. 5. The Hamiltonian and noise parameters taken in 
this study reproduce the experimentally-observed absorption 
spectrum. 


In this study we have employed a model of the R. 
rubrum LH1 where 32 BChl molecules, bound to 16 a and 
fd polypeptides as a/ 3 (BChl )2 subunits, are equidistantly 
arranged along a ring of radius r = 47 A. Dipole moments 
are taken to be tangential to, and in the plane of, the ring. 
The sub-nanometer distance between neighbouring chro- 
mophores implies that the nearest-neighbour couplings 
depend on the geometry of the electronic wave-function 
of each chromophore and cannot be inferred directly from 
a dipolar interaction in the point-dipole approximation. 
Couplings on the ring can be fitted using fluorescence 
anisotropy measurements, resulting in nearest-neighbor 
couplings Q\ = 600 cm ' 1 and Q 2 = 377 cm -1 for the 
intra- and inter-dimer couplings, respectively^. As dis¬ 
cussed in the text, homogeneous broadening is taken into 
account using either a Lindbladian master equation with 
a pure clephasing rate of 7 d = 106 cm ' 1 or the equiv¬ 
alent protocol of dressing excitonic stick-spectra with a 
Lorentzian of full-width half-maximum 2 x 7 ^. These 
parameters reproduce the experimentally-measured ab¬ 
sorption spectrum of isolated LH1, see Figjbj Coupling 
between dipoles on different rings, separated by a = 120 
A, is calculated using the dipole-dipole interaction 

f ( dm ■ d n ^ i^dn ' Ar nm ^ ^ 'd m * ^r nr 2j 
J ”’ m “ 4 ttk |Af„, m | 5 j 

(Al) 

where Ar njm = f n — r m and «=2 is the relative permit¬ 
tivity. The LH1 dipoles are taken to be of magnitude 
y/2A x 6.3 D = 9.8 D to take into account the observed 
hyperchromism^. 


Appendix B: Measures of extended delocalization 


Excitonic properties from purple bacteria membranes 
can be obtained from the Coulomb exchange Hamiltonian 

2NR 

H = ^2 ( e m + $e m ) I m) (m\ 

m 

2NR 

+ ^2 J mn{ \n) ( m\ + |to) (n|) (Bl) 

m^n 

where J nm denotes the interaction strength between the 
induced transition-dipoles of pigments n and m, corre¬ 
sponding to electronic excited states |n) and | m) with 
energies e n and e m . Note that each ring contains 2 N = 32 
pigments and we examine a line of rings (Q=l), as in the 
main text. The wavefuntions take a remarkable simple 
form in the absence of protein inhomogeneities {5 tm = 0), 
as given in eq (5) of the main text. 

Explicitly, the bright states of relevance for the optical 
response are 

I ±,k x = k v = 1,Q = 1) oc ^]sin l±> r ) • 

(B2) 


These states present populations in each ring propor¬ 
tional to sin 2 ( 7 rr/(R + 1)) which is largest for the rings 
lying in the middle of the linear chain. Note that in 
the abscence of coherent interaction among rings, all 
states should have populations concentrated in single 
rings. As a consequence, on average, a density operator 
will present a statistical mixture with populations being 
equally shared among all rings. 

Absorption is a linear functional in the density op¬ 
erator. Hence, we concentrate on the density operator 
p = | a) (a | which results from averaging the eigenstates 
of eq (B 1 ) under realistic conditions, i.e. with appropri¬ 
ate magnitude for protein-induced fluctuations y/(6%) = 
325 cm -1 . Our choice of magnitude of these fluctuations 
is justified by a detailed analysis of absorption spectra of 
LH1 in R. rubrum, see Fig. [5j 

The average population of each ring for the subset of 
states |a) with the largest light-induced dipole strength, 
is presented in Figure [ 6 ] Here, a non-uniform population 
in the linear configuration of rings can be seen, where 
the complexes in the middle of the chain present a greater 
population compared to the rings at the chain edges. The 
pattern of ring populations displayed in the presence of 
inhomogeneities resembles the noiseless case from states 
|±, k x = k y = 1, Q = 1), which illustrates the formation 
of excitons fulfilling the full array boundary conditions in 
the presence of quasi-static fluctuations. The resilience 
of the excitons to the inhomogeneous disorder illustrates 
that the disorder produce noticable smearing out of ex¬ 
citonic delocalization, but does not completely erase its 
features. Another indicator of excitonic coherence among 
neighbouring rings is also presented in Figure [ 6 ] namely 
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FIG. 6. The linear array of rings in top presents ring populations ^2i nfLr \{n\ p\n)\ an( l total inter-ring coherences 
X/nei-! mer 2 I ( n l P I m ) I shown in black and red bars, as a result from averages of 50000 stochastic realisations of inhomoge¬ 
neous noise. The result of the the noise-less model for populations is presented in blue connected crosses. 


the total coherence p ri , r2 = J2nen,mer 2 I ( n l P \ m ) I among 
all chromophores that belong to neighbouring rings n and 
r 2 - This total inter-ring coherence, ussually assumed to 
be negligible for inter-ring exciton transfer, is manifestly 
non-zero in Figure [6] and is the cause for the absorption 
anisotropy that produces the finite linear dichroism ( LD ) 
discussed in the main text. 

To characterise the excitons formed upon photon ab¬ 
sorption, useful quantities are the conventional inverse 
participation ratio 


IPRc = 



(B3) 


and a generalisation which we denote the ring inverse 
participation ratio 


IPR r = 




(B4) 


tend well beyond the current paradigm of single ring ex¬ 
citons (lying flat with IPRr = 1), sufficient for the de¬ 
scription of processes with a longer built-in timescale. 

The dependence of excitons length with the size of the 
array is made conspicuous in Figure [7](b) where the av¬ 
erage IPRr is shown as a function of the linear (Q = 1) 
array size for the two brightest states |a) of each realisa¬ 
tion. Here, it can be appreciated an increase in exciton 
size with the array extension in the same qualitative man¬ 
ner as the homogeneous wavefunctions of eq (B2 ) (also 
shown in the same figure). 

Based on the observations that the wavefunction and 
the excitonic length from averages of inhomogeneous 
noise resemble qualitatively the features that arise from 
a noiseless treatment, it can be safely assumed that the 
inhomogeneous broadening is included, the noiseless fea¬ 
tures are not completely smeared out. It becomes natu¬ 
ral therefore, to assert the properties of the system based 
on the noiseless characterization presented through equa¬ 
tions (5)-(7) in the main text. 


where the inner and outer sums are performed over BChls 
that belong to a specific ring r and over all rings in the 
array, respectively. 

As defined, IPRc ranges from 1 to the total num¬ 
ber of pigments in the array and measures how many 
chromophores participate in a given exciton, while IPRr 
ranges from 1 for any exciton confined to a single ring, 
to IPRr = R for pure states that are evenly delocalised 
over the entire array of R rings. For IPRr > 1 or 
IPRc > 32, these two quantities demonstrate unam¬ 
biguously the presence of bright excitons over domains 
greater than a single ring (see Figure |7]a)). A general 
trend shows larger exciton lengths for arrays with higher 
number of rings. Figure UN highlights realisations ex¬ 
tending over more than 2 rings ( IPRr > 2), which ex¬ 


Appendix C: Shift of the unpolarised absorption spectrum: 
dependence of Vi,i on the geometry of the ring dipoles 

The resonant interaction betwee rings in general as¬ 
semblies can result in a shift of the absorption spectra, 
addressed by equation 6 in the main text and numerically 
supported by Figure 2(b) and 2(d) in the main text. The 
magnitude of this shift, in comparison to the isolated ring 
situation, depends on the structural details of the rings. 
As it has been stated in section [A] of this SI, we have 
used a structure where the Q y transition dipoles are tan¬ 
gent to the circumference of the ring. However if this 
angle, Ay, is increased, the magnitude of increases, 
as presented in Figure [8j Accordingly, the determination 
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FIG. 7. (a) Scatter plot (points) of the IPRc and IPRr from 
eqs.( B3|B4| | for two (red) and five (blue) rings, for states with 
the greatest dipole moment in realisations of the quasi-static 
fluctuations in eq |B1| States that delocalize beyond 2 rings 
are highlighted. In (b) is presented the average IPRr for the 
brightest state in each noise realisation (crosses) and for the 
noiseless wavefunction in eq (B2| (circles), as a function of 
the length of the array. Symbols connected by straight lines 
to guide the eye. 


of the shift in the spectra reflects specific details of the 
rings, which escape measures performed on single rings, 
but are otherwise accessible through the observation of 
assemblies. 


Appendix D: Inhomogeneous broadening. 

The spectral lineshape and width are determined by 
both homogeneous and inhomogeneous contributions. 
Two effects are manifested when inhomogeneities are in¬ 
troduced in excitonic transitions, namely, lineshape mod¬ 
ifications and energy shifts. Figure[9](a) shows that over- 



FIG. 8. Coupling Vi,i as a function of the angle Ay between 
the Q y transition dipole and the ring tangent. 


lapping spectra are obtained with different combinations 
of homogeneous and inhomogneous broadening, being the 
most conspicuous signature of their slightly different line- 
shapes, the difference in the low energy tails. For the case 
of smaller inhomogeneous broadening, a Lorentzian-like 
lineshape appears with slower rise as compared to the re¬ 
sult with higher inhomogeneity with a Gaussian-like line- 
shape. For values of inhomogeneous disorder bounded by 
the full spectral width, the energy splitting between ab¬ 
sorption spectra along orthogonal polarizations is greater 
with the inhomogeneous noise Figure [9] (b). This phe¬ 
nomenon can be understood by the excitonic repulsion 
which results in slightly larger splittings for greater lev¬ 
els of nois^H. Accordingly, the LD changes due to the 
inhomogeneous noise (see inset of Figure [9] (b)), with a 
very counter-intuitive result, namely, greater inhomogen- 
ities produce greater LD amplitude. Such increase can 
only be traced back with full numerical simulations, as it 
requires to account for all the dark levels in order to quan¬ 
titatively assert the shift calculated in the inset of Figure 
[9] (b). Note however that within the inhomogenities ex¬ 
pected in these systems, the energy difference among or¬ 
thogonal absorption spectra does not produce dramatic 
changes in the LD , given that we have been careful to 
reproduce in all cases the absorption profile shown in Fig¬ 
ure [9] (a). The major change in the LD will be due to 
the full broadening of the optical transition transition - 
labeled in Figure [9] (a) as I - as it becomes conspicuous 
by comparison of the LD from homogeneously broadened 
spectra with a width of 2x106 cm -1 in Figure 1 of the 
main text, and the LD of full spectrum, here presented. 
In the former case, the amplitude of the LD can reach 
about 5% of the absorption spectrum, while in the lat¬ 
ter case the LD amplitude is reduced quite generally to 
2 — 3%, as Figure [9] (b) illustrates. It can be concluded 
from this analysis, that the LD changes only slightly de¬ 
pending on the inhomogeneous component in the spectra 
with a robustness only vulnerable to the total width T of 
the absorption spectrum. 

The noise model presented above accounts for 
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FIG. 9. Effects of energetic and structural disorder on LD. In (a), the overall absorption spectra takes on more Gaus¬ 
sian lineshape, manifested from steeper tails of the spectra as inhomogeneous disorder is increased. The splitting between 
orthogonally-polarized manifolds, A E, in (inset b), has only a slight increase with inhomogeneous noise, which results in an 
LD signal (b) robust to energetic disorder (ainhom = 175 - 325 cnG 1 ). The LD is also robust to variations of the inter-ring 
distance, as shown in (c). In (d) the LD ensemble averaged over randomly-oriented oval structures in linear arrays with 1, 2 
and 3 rings. 


environmentally-induced inhomogeneities in the pigment 
energies. However, it does not account for explicit pertur¬ 
bations of the chromatophore vesicle, such as variations 
in the inter-complex distance or for structural perturba¬ 
tions, such as elliptical deformations of the rings. 

Examination of paracrystalline 67 domains of LH2 in 
Rsp. photo metricum and close-packed orthorhombic 
crystals2§--®l in R. rubrum imply variations of only a 
few Angstroms in the cente r-to-c enter distances among 
neighbouring ring structure :] 67 ^ 70 ! Figure |9](c) shows the 
effect on the LD when the inter-complex distance suffers 
stochastic variations which are much greater than the 
variations observed from pair correlation functions calcu¬ 
lated from atomic force microscopy (AFM) experiments 
(a/ ( 5r 2 ) = 0,5,10,15 A as compared to yj ( Sr 2 ) ~ 1 
A observed iiP). This figure illustrates that the con¬ 
trast in the LD actually increases with greater center- 
to-center distance variations. This surprising result can 
be understood from the functional form of the Coulomb 
interaction. At a center-to-center distance of a few tens 
of nanomenters, the Coulomb interaction scales as 1/r 3 
which means that bringing rings together by Sr increases 
the coupling more than it is decreased by separating them 


by Sr, hence, the LD is enhanced in average. A fur¬ 
ther signature is apparent for large variations reflected 
by the red shift of the LD. This signature can be traced 
back to the amplitude of absoprtion along x and y axes, 
which for strong variation of the center-to-center distance 
makes the former absorption greater than the latter, and 
hence produces a LD which is slightly asymmetric. Note, 
however, that the few Angstrom variations observed in 
experiments are too small to appreciably affect the LD 
signal by its magnitude or to enhance appreciably the 
absorption along one of the polarization direction with 
respect to the orthogonal, by means of inter-ring distance 
inhomogeneities. 

In general, LH complexes in purple bacteria express 
geometries which may differ from the circular ring of 
R. rubrum in our model. Another conformation of the 
LH1 in R. rubrum has been observed to be significantly 
ellipsoidal. In addition, the LH2 of Rubrivax gelati- 
nosus, obtained by AFXp3, and the latest LH1 X-ray 
structure of Tch. tepidum both exhibit a degree of 
ellipticity^. It is not known if the long axes of ellipti¬ 
cal rings may be aligned or randomly oriented along the 
membranes where they are synthesized. A study of el- 
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liptical structures is therefore in order as it might clarify 
the scope of the LD to witness extended delocalization. 

Oval structures with an apparently small eccentricity 

yjl — (|) 2 = 0.32, where a and b are the minor and 
major axes, respectively, produce energy shifts among 
k = ±1 states of ~ 100 cm -1 in a single rin^H, which is 
10 times larger than the shifts due to inter-ring excitonic 
coupling, presented in the inset of Fig. 5 (b). In these 
experiments the excitation polarisation was fixed with re¬ 
spect to the major axis of the single elliptical ring. In an 
ensemble, there is no reason to believe that the orienta¬ 
tions of major axes will or will not correlate. In order to 
isolate the effect from elliptical deformation, we momen¬ 
tarily set the pigment energy inhomogeneities to zero. 
Now, if an ensemble presents randomly-oriented major 
axes and the rings do not couple coherently, then the 
finite LD signals, stemming from the elliptical geome¬ 
tries, cancel out and lead to a vanishing LD (see Fig. [9] 
(d)). In order to account for the elliptical deformations 
we use an interaction between nearest neighbours with a 
functional form that only depends on their relative dis- 

tance J n ,m — Q 1.2 where ro ~ 9.2 A is the distance 
among pigments in the circular geometry, and Qi ,2 are 
the nearest-neighbour couplings. In the same fashion, if 
an ensemble of randomly-oriented coupled rings is consid¬ 
ered, then the LD arising from the elliptical perturbation 
vanishes and only the nonzero signal arising from the co¬ 
herent interaction remains. Figure 4(e) shows that the 
nonzero LD for coherently-coupled linear arrays is robust 
to the ensemble-averaged elliptical perturbation, as it re¬ 
sults in a maximum value of 5 % for two rings, which, as 
expected, is similar to the value, with just homogeneous 
noise, presented in Fig. 1 of the main text. 

It was observed that the macromolecular assembly in 
purple bacteria vesicles can indeed present long-range 
organisation 73 . Therefore, the likelihood of having affin¬ 
ity domains in each ring that set specific connecting 
residues to neighbouring units and result in an overall 
specific orientation of all major axes, is not out of the 
question. If the macromolecular assembly has a preferred 
orientation of rings-ellipse axes, a finite LD will be mea¬ 
sured even in the absence of appreciable inter-ring exci¬ 
tonic delocalization. 

In general, if the geometrical contribution to the LD 
is independent of the size of the assembly, then a sim¬ 
ple subtraction among LD signals, A LD, arising from 
different assemblies lengths will therefore represent the 
residual contribution resulting from the extended exci¬ 
tonic delocalisation among such assemblies. 


Appendix E: General 2D geometries 

The description in the main text of excitonic delocal¬ 
isation over extended domains made use of square unit¬ 
cell lattices. It was concluded that arrays with a large 
aspect ratio and small width are desirable in order to 
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FIG. 10. LD spectra for the triangular lattice configuration 
and increasing array size. Insets show the lattice configura¬ 
tion studied and stick polarized specta. A small asymmetric 
array, as in (a), exhibits a larger LD than its double (b) due 
to a decreased splitting of polarized states. This can be qual¬ 
itatively understood from the results for the square unit-cell, 
eq (7) of the main text. 


obtain an polarized optical response which encodes the 
symmetries of the extended excitons. Here we generalise 
these statements to arrays with other unit cells. Specif¬ 
ically, we analyse a triangular lattice and, to prove the 
generality of our statemements, we also study a square 
unit cell with non-nearest-neighbour interactions. 

The possibility to use the natural assembly of these 
complexes to probe inter-ring excitonic delocalization is 
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investigated by numerical diagonalization of the Hamilto¬ 
nian in the manifold fc|. | k'\ = 1 for a triangular unit-cell 
lattice (Figure [l0|. The triangular unit cell approximates 
the natural “quasi-hexagonal” aggregation state observed 
for LH1 complexes in photosynthetic membranes of R. 
rubrurr^. Triangular para-crystaline domains of LH2 
complexes also form in Rsp. photometricum at low light 
intensity conditional. Figure flOfa), shows that the dipole 
strength in the triangular lattice is concentrated at the 
band edges, analogous to the prediction made by eqs (3- 
5) in the main text. As shown in this figure, states po¬ 
larized along and perpendicular to the long-axis of the 
array are split in energy, leading to a finite LD once the 
stick spectra is dressed with appropriate line-shape func¬ 
tions. Figure |To| (b) highlights another important feature. 
Here, even though an asymmetric array is being consid¬ 
ered with the same aspect ratio, the LD signal is, in fact, 
smaller ( LD max ss 0.6 % in (b) as opposed to LD max ss 
0.7 % in (a)). This can be understood from eq (7) in the 


main text, which shows that the energy splitting between 
bright states decreases for arrays of the same aspect ra¬ 
tio but increasing size. Hence, independent of unit-cell 
configurations, only arrays with a width of a few rings 
will be suitable to probe extended excitonic delocaliza¬ 
tion through LD. 

The same conclusion can be drawn based on a square 
unit-cell with non-diagonal interactions, as shown in 
Figure El The additional interactions produce addi¬ 
tional splittings, and arrays with appreciable aspect ratio 
( Q/R 1 or Q/R <C 1) still present a dipole moment 
redistribution that leads to polarization anisotropy. If 
the array is symmetric, the polarized states are not split 
in energy and the array shows zero LD , which can be 
seen from eq (7) in the main text. These results show 
that, independent of the unit-cell geometry, it is possi¬ 
ble to witness LD in general arrays of small width and 
appreciable aspect ratio. 
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FIG. 11. LD spectra for the square lattice configuration and 
non-nearest-neighbor interactions. Insets show the lattice 
configuration studied and stick polarized specta. The strip 
of dimers in (a) reflects the native organization of LHls in 
membranes containing LH1 and LH2, and exhibits a non-zero 
LD. Symmetric arrays (Q = R), as in (b), show zero LD since 
the polarized states are degenerate in energy. These results 
are qualitatively similar to those for the square unit-cell with 
nearest-neighbor interactions, eq (7) of the main text. 

























